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We study a Bose Josephson junction realized with a double-well potential. We propose a strategy 
to observe the coherent superpositions of phase states occurring during the time evolution after a 
sudden rise of the barrier separating the two wells. We show that their phase content can be obtained 
by the full counting statistics of the spin-boson operators characterizing the junction, which could 
be mapped out by repeated measurements of the population imbalance after rotation of the state. 
This measurement can distinguish between coherent superpositions and incoherent mixtures, and 
can be used for a two-dimensional, tomographic reconstruction of the phase content of the state. 
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I. INTRODUCTION 

Superpositions of quantum states are at the heart of 
quantum physics. They are important for precision tests 
and arc at the basis of quantum information and commu- 
nication. The fate of these superpositions in the presence 
of interactions with the environment has been largely 
studied to understand the cross-over from a quantum to a 
classical description Although fragile and subject to 
decoherence, superposition states have been experimen- 
tally observed with trapped ions and cavity photons 
1- 

Ultracold atomic gases are another system potentially 
suitable for realizing macroscopic superpositions. We are 
interested here in the study of superpositions of (angular- 
momentum) coherent states 0], which could be realized 
in a Bose Josephson junction made by a double- well trap. 
Such superposition states are predicted to form [f| during 
time evolution after a sudden rise of the barrier separat- 
ing the two wells, starting from a coherent state at the 
initial time. An initially coherent state corresponds to 
the ground state in the regime where the tunneling en- 
ergy dominates the repulsion energy, and is commonly 
realized in current experiments on Bose Josephson junc- 
tions [6(. The subsequent dynamical evolution after the 
rise of the barrier is driven by interactions only, and at 
specific times the system undergoes the formation of sev- 
eral multicomponent superpositions of phase states, fi- 
nally coming back to the initial state after a revival time. 
In a pioneering experiment 0] on ultracold bosons in op- 
tical lattices, revivals have been indeed observed, but it 
was not possible to directly access the macroscopic super- 
positions which were formed during intermediate stages 
of the evolution. 

Multicomponent superpositions of phase states are 
characterized by a sequence of equal and equidistant 
peaks along the equator 8 = tt/2 of the Husimi phase 
distribution Q{0,<f>) = (6<t>\p\6<t>) @, H, with \6<f>) an 
angular-momentum coherent state and p the density ma- 



trix. Each peak corresponds to the contribution of a 
given coherent state belonging to the coherent superpo- 
sition. An important issue to address is how to reveal 
the effect of decoherence on quantum superpositions of 
macroscopically distinguishable states. We proceed by 
analogy with the macroscopic superpositions of phase 
states studied in the context of photons in cavities [Io| . 
We expect that generic sources of noise would tend first 
to destroy the coherence between the components of the 
superposition, by projecting it onto an incoherent mix- 
ture of coherent states. On a longer time scale each com- 
ponent is expected to relax, yielding a total spreading of 
the phase profile. 

In this work we address the question on how to identify 
the quantum macroscopic superposition of phase states 
and how to distinguish them from mixtures. Our ap- 
proach is substantially different from the one of Ref. [11| 
which is devoted to map out the Husimi distribution. 
Although the Husimi phase distribution is in one-to-one 
correspondence with the full density matrix in prac- 
tice it is almost insensitive to the difference between a 
coherent superposition of phase states and the corre- 
sponding incoherent mixture. This is because the Husimi 
distribution is the diagonal of the density matrix repre- 
sented over coherent states. Figure Q] (top panel) shows 
the Husimi distribution for a three-component superpo- 
sition of phase states and for the incoherent mixture, the 
tiny difference between the two being illustrated in the 
inset. 

The paper is organized as follows. After introducing 
the model in Sec. II, in Sec. Ill we propose to reconstruct 
the phase distribution by mapping out all the higher or- 
der correlations, ie the full counting statistics of the spin- 
boson operators J x and J y associated to the quantum 
model of the junction. In Sec. IV we show that it is 
possible to make a two-dimensional reconstruction of a 
Wigner-like distribution function. We conclude in Sec. 
V after discussing experimental issues. 
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FIG. 1: (Color online) Top panel: Dimensionless Husimi dis- 
tribution Q(9 = 7r/2,0) for a three-component superposition 
of phase states (solid line) and for the corresponding incoher- 
ent mixture (dashed line), as a function of the phase <j)/2n for 
N=20 particles. The inset shows a zoom of the same function 
around <j> — 2tt/3, illustrating the difference between the su- 
perposition state and the incoherent mixture. Bottom panel: 
section of the Bloch sphere on the equatorial plane 9 = it/2, 
parametrized by the angle <j>- The dots indicate schematically 
the phase of the three coherent states which give rise to the 
superposition, and correspond to the maxima in the Husimi 
distribution illustrated in the top panel. The vector r defines 
a generic direction of the spin-boson operator for which the 
probability distribution is considered (see text). Our conven- 
tion for the x,y axes is also indicated. 



II. TWO-MODE MODEL FOR THE BOSE 
JOSEPHSON JUNCTION 



We describe the Bose Josephson junction in the quan- 
tum regime by using the two-mode approximation. This 
model applies in the limit where the chemical potential 
of the system is smaller than the barrier between the two 
wells and corresponds to neglecting higher order single- 
particle levels of the double-well potential. Hence we 
adopt the following two-mode Bosc-Hubbard Hamilto- 



a A ft t * ^ 1 t i 

H = Eia[ai + E2a\ai + — a[a[aidi 
+ ^-a\a\a,2a2 - K(a\a\ + a\a 2 ), 



(1) 



where a,i,a} with i = 1,2 are bosonic field operators for 
the bosons in each well, satisfying the usual commuta- 
tion relations [a^, a]] = 5ij; Ei are the energies of the two 
wells, Ui > are the boson-boson repulsive interactions 
and K is the tunnel matrix element, ie the Rabi oscilla- 
tion energy in the case of a non-interacting model. We 
then perform the transformation onto spin-boson vari- 
ables, where J x = (a\a,2 + a,2CLi)/2 is the inter-well tun- 
neling operator, J y = — i(a\a,2 — «2ai)/2 is the current 
operator, and J z = (a[ai — a\a,2)/2 = n is the popu- 
lation imbalance operator. The Hamiltonian is readily 
rewritten as 



H = U S {J Z -n ) 2 /2-2KJ x 



(2) 



with U s = (Ui + U2) and no is related to the imbalance 
between the wells according to n = —\E\ + U\{N — 
l)/2 - E 2 - U 2 {N - l)/2]/U s 0. In the quasi-classical 
limit KN ^> C/ s , the angular-momentum coherent states 

1/2 



Mo) = E 



N/2 



with 



N 

--n/2 y m + N / 2 J (i+M 2 )"/ 2 

a = tan(# /2) exp(— icj) ) are useful to describe the state 
of the system. In particular, the ground state corre- 
sponds to the coherent state with 6q = 7r/2 and (f>o = 0. 

In the following we consider the time evolution of the 
Hamiltonian ([2|) after a sudden rise of the barrier, namely 
we set K = 0, starting from the nonintcracting ground 
state \9q = 7r/2,0o = 0) as the initial state. Such a dy- 
namical evolution occurs entirely in the equatorial plane 
of the Bloch sphere as J z is a constant of the motion. 
We focus on the specific time T q = 2tt /(U s q) where a 
coherent superposition of q coherent states is formed, ac- 
cording to 



q-l 
fe=0 



(3) 



with uq = ie 2 « Yftr- 



,1,, 



-1) 



we have assumed q odd. Our aim is to reveal the phase 
content of such a state. For this purpose we consider 
the operators J x and J y which have an action close to 
that of the relative-phase operators cos <j) and sin </>. How- 
ever, the latter act incorrectly on the boundary states 
\±N/2). The average over a coherent state characterized 
by a relative phase 4>q and population imbalance n gives 
(a\J x \a) = ncos4>o and (a\J y \a) = Tisin^o. While the 
momentum distribution accesses only the first-moment 
averages (cos^>), (sin^>) [13], in order to access the en- 
tire phase distribution we need to consider all the higher 
correlations. 
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III. FULL COUNTING STATISTICS OF THE 
ANGULAR-MOMENTUM OPERATORS 



The probability distribution of the eigenvalues of an 
observable contains the information equivalent to the 
knowledge of all the moments of the distribution itself. In 
turn, the latter are known once the generating function 
is known, according to a full counting statistic approach 

EH- 

Let us consider a spin-boson operator J r in the di- 
rection specified by the two-dimensional vector r = 
(sin<^>, — cos0) in the (x, y) plane. The generating func- 
tion for the probability distribution of the eigenvalues r 
that correspond to the eigenstates \r), J r \r) = r\r), is 
given by 



for the incoherent mixture and 



-irfj r 



(4) 



Here, J r = J x sin0 — J y cos0 and (...) indicates the quan- 
tum average over the state of the system. Considering for 
simplicity a pure state \ip) and expanding in terms of the 
eigenstates \r) of J r we have 



M 7 ?) =X/ 



-%r\r I 



<V>|r)| 2 



(5) 



where the eigenvalue r takes integer values in the inter- 
val [—N/2, N/2] and we have assumed N to be even. 
We shall denote P<p{r) = \{tp\r)\ 2 the corresponding- 
probability distribution (for a generic density matrix p, 
P^{r) = Tr{p\r)(r\)). Its shape reflects the phase content 
of the state projected along the direction specified by the 
vector f (see the bottom panel of FigQ]for a sketch). As 
we will detail in Sec.V, this could be accessed experimen- 
tally. 

In order to determine the probability distribution 
P<j>{r), we evaluate the generating function Eq.Q ana- 
lytically both for coherent superpositions and mixtures. 
Specifically, we focus on the three-component superpo- 
sition of phase states and on the corresponding incoher- 
ent mixture of the same three equally-weighted phase 
states. We use the definition of angular-momentum co- 
herent states and the disentangling formula [H 



hf h (v) = hf xt ( V ) 

9-1 - ~* 
CfcC fc , 



(8) 
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for the coherent superposition. The probability distri- 
bution (r) is readily obtained by Fourier transforming 
the above expressions. 

By considering the projections along the x and y di- 
rections, respectively, wc obtain the full counting statis- 
tics of the operators J x and J y . This is illustrated in 
FiglUfor the three-component superposition |-0 3 ) of phase 
states as well as for the corresponding mixture. The 
distribution is peaked around the semiclassical values 
for (j x ) = (JV/2)cos(±7r/3), {N/2) cos(tt) and {J y ) = 
(JV/2)sin(±7r/3), (JV/2) sin(?r). We find that the distri- 
bution P-n/2^) displays a noticeable difference between 
the mixture and the coherent superposition: the latter 
displays oscillations which are absent in the former. The 
presence of fringes in the distribution of the eigenvalues 
of angular momentum operators for superposition states 
was also noticed in the context of the dynamics of the 



quantum non-linear rotator by Sanders [lj]. The func- 
tion P-k{t) instead does not display fringes for the three- 
component superposition because its components do not 
overlap when projected along the y-direction (see Fig[TJ 
bottom panel); as a result no interference effect takes 
place in this case. 

The result for the three-component state extends to 
higher-component superpositions; the two-component 
one instead cannot be distinguished from the correspond- 
ing incoherent mixture by this method, due to the specific 
form of its state components. The full counting statistics 
of the operator J z could also be defined, but does not 
yield any useful information about the considered super- 
positions of phase states as it coincides with the binomial 
distribution P^oM = ^r(il^_ r ) of the initial coherent 
state. 



e -i>lJr =e -rV_ e -log(l+|r| )J, e rJ +) (g) 



with r = tan(r?/2)e % ^ . We thus obtain the final results 
for the generating functions: 



9-1 

h^{ V ) = \u \ 2 J2 
k=0 



(7) 
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IV. TWO-DIMENSIONAL RECONSTRUCTION 
OF THE PHASE DISTRIBUTION 

A two-dimensional (2D) tomographic reconstruction of 
the phase content of a state in the {x, y) plane is possi- 
ble, using the concept of the Radon transform [15^ . This 
has been also used for coherent superpositions with cav- 
ity photons The 2D distribution function f{x,y) is 
obtained once all the one-dimensional projections P<p{r) 
are known as a function of <j) in the interval [0, 2ir]. One 
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FIG. 2: (Color online) Eigenvalue distribution P$(r) corre- 
sponding to J x (<f> — 7r/2) and J y (4> = it) for the three- 
component coherent superposition (solid lines) as well as for 
the incoherent mixture of the same phase states (dashed lines) 
with N = 20. The vertical lines correspond to the semiclas- 
sical values for (j x ) and (J y ) for the coherent states entering 
the superposition. 



should invert the expression 

/+oo /- + OC 

dx I dyf(x, y)S(r — x sin (j> + y cos (/>), 
-oo J — oo 

(9) 

which, using the definition of the generating function in 
Eq. |5j), yields as a final result 



f(x,y) 



(27T) 2 



T}dr) I d(j>h^(ri)e 
Jo 



irj(x sin (fi — y cos c; 



(10) 

This is a quasi-probability distribution for the non- 
commuting operators J x and J y . It could be regarded as 
the two-dimensional projection on the equatorial plane 
of the Bloch sphere of a SU(2) Wigner function in the 
spirit of Ref . [17[ . It is readily verified that the marginal 
probability distribution obtained by integration of f(x, y) 
along one of the two directions y or x yields the one- 
dimensional distributions P(p(r) with <fi = 7r/2 and 7r, 
respectively. The function f(x,y) is not in one-to-one 
correspondence with the state of the system. However, 
for the specific superpositions of phase states which we 



consider here it yields the main information about the 
phase structure of the state. Furthermore, in the case of 
a coherent superposition of phase states it is markedly 
different from the corresponding incoherent mixture. 

Figure [3] illustrates the 2D quasi-probability distribu- 
tion f{x,y) for a three component coherent superposi- 
tion. It shows three pronounced maxima in correspon- 
dence with the three coherent states giving rise to the 
macroscopic superposition. It also displays oscillations 
between the maxima, due to interferences between the 
components. The 2D quasi-probability function evalu- 
ated for the corresponding incoherent mixture also ex- 
hibits the main peaks but the fringes are strongly sup- 
pressed (see FigJ3]bottom panel), the small remaining os- 
cillations being intrinsically due to the definition of the 
function f(x,y) as a Fourier transform in angular vari- 
ables. 



V. EXPERIMENTAL ISSUES AND 
CONCLUDING REMARKS 

For each choice of the angle <j> the probability distri- 
bution P^(r) can be experimentally accessed by repeated 
measurements of the corresponding angular momentum 
operator J r . Indeed, since the eigenstates of J r form an 
orthonormal basis, each superposition state decomposes 

as \%) = J2ri-N/2 c r\ r ) witn c r = (»#«)• Then, accord- 
ing to the postulates of quantum mechanics, after a (pro- 
jective) measurement of J r the state jumps to the state 



with probability PS° (r) 



r.912 



, and the correspond- 



ing outcome of the measurement is r. The full distribu- 
tion P<f,{r) is obtained by repeating this procedure many 
times, each time preparing the system in the same initial 
state [3]. The measurement of J r for a generic angle <f> 
can be achieved by measuring the population imbalance 
J z between the two wells - a variable typically accessed 
in experiments @ - after proper rotations of the state 
over the Bloch sphere [2(|. More precisely, measuring 
J r on the quantum state l^) is achieved by measuring 
J z on the rotated state e I7rJa= / 2 e^ Ji I?/)). Each rotation 
can be implemented by a time evolution of the system 
with Hamiltonian j2]) once the interaction strength U s is 
tuned to zero (eg by a Feshbach resonance). The rotation 
direction and angle can be chosen by tuning the height 
of the barrier, the ener gy b ias between the wells and the 
time-evolution interval [2l|. 

In summary, using the concepts of full counting 
statistics we propose a method to access to the full 
phase content of a macroscopic superposition of phase 
states, as well as to reconstruct tomographically its two- 
dimensional quasi-distribution. This method is capable 
to distinguish between such states and an incoherent mix- 
ture as the latter does not display fringes in the probabil- 
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FIG. 3: (Color online) Top panel: Dimensionless 2D quasi- 
probability distribution f(x,y) in the (x,y) plane (dimen- 
sionless) for a three-component coherent superposition with 
N = 20 particles. Bottom panel: section of the same quasi- 
probability distribution (solid line) in the same units, in 
the direction y = 0. The dashed line represents the quasi- 
probability function f(x,0) for the corresponding incoherent 
mixture of the same three coherent states. 



ity distribution. The superposition of phase states could 
be realized with current experiments on Bosc Joscphson 
junctions. 

Feasibility of the current proposal requires control of 
the atom number on the experiment; atom losses are pre- 
dicted to destroy the superposition and seem to present 
the major experimental challenge. 

During the preparation of the present manuscript we 
have became aware of similar results for the 'NOON' 
state by Haigh et al. [22| . 
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